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ABSTRACT 

We report the detection of shock-excited far-infrared emission of H 2 O, OH, and CO from 
the supernova remnant 3C 391, using the ISO Long-Wavelength Spectrometer. This is the 
first detection of thermal H 2 O and OH emission from a supernova remnant. For two other 
remnants, W 28 and W 44, CO emission was detected but OH was only detected in absorption. 

The observed H 2 O and OH emission lines arise from levels within ^ 400 K of the ground 
state, consistent with collisional excitation in warm, dense gas created after the passage of the 
shock front through the dense clumps in the pre-shock cloud. The post-shock gas we observe 
has a density ^ 2 x 10® cm“® and temperature 100-1000 K, and the relative abundances of 
C0:0H:H20 in the emitting region are 100:1:7 for a temperature of 200 K. The presence of a 
significant column of warm H 2 O suggests that the chemistry has been significantly changed by 
the shock. The existence of significant column densities of both OH and H 2 O, which is at odds 
with models for non-dissociative shocks into dense gas, could be due to photodissociation of 
H 2 O or a mix of fast and slow shocks through regions with different pre-shock density. 


Subject headings: supernova remnants - Infrared: interstellar: lines - ISM: molecules - ISM: 
individual (3C 391) 


1. Introduction 


Long suspected to be an important component of the interstellar medium, water has proven difficult 
to observe because of the dominating absorption by water vapor in the Earth’s atmosphere. Water masers 
at radio frequencies are bright enough to penetrate the atmosphere, but they are rare and their physical 
conditions may be exceptional. Now that we have flown a sensitive infrared spectrometer in space, aboard 
the Infrared Space Observatory (Kessler et al. 1996), transitions among energy levels ~ 100-1000 K 


^Based on observations with ISO, an ESA project with instruments funded by ESA Member States (especially the PI 
countries: France, Germany, the Netherlands and the United Kingdom) with the participation of ISAS and NASA. 
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above the ground state of H 2 O are observable. Early results from ISO indicate that H 2 O is a significant 


constituent in a cloud near the galactic center ( 

Cernicharo et al. 1997 

) and that H 2 O is as abundant as CO 

in shocked regions in Orion ( 

Harwit et al. 1998 

)• 



The observations reported here are part of a study of the infrared emission from supernova remnants 
interacting with molecular clouds. Our targets were selected from a sample of supernova remnants with 


previous evidence for interaction with nearby molecular clouds based on X-ray and radio morphology (Rhc 
fc Petre 199g ), millimeter-wave molecular line observations ( |Wilner et al. 199E , Wootten 1977 , Wootten| 
Reach and Rho 1998 ), and OH 1720 MHz observations ( Frail et al. 1996 , ^reen et al. 1997 ). The 


1981 


presence of 1720 MHz maser emission without bright main-line OH maser emission suggests that they are 
collisionally excited ( Elitzur 1976| ), and their location in bright radio remnants indicates that they are 
produced in the dense regions just behind molecular shocks. The first results of our ISO observations 
showed that the [O I] 63 /rm line, which is expected to be one of the brightest cooling lines for post-shock 
gas for a wide range of gas densities and shock velocities, is very bright at the locations of the OH masers 
( Reach and Rho 1996 ; Paper I). In a parallel observational project, we also found shock-accelerated CS, 
CO, and HCO+ molecules in 3C 391 (Reach and Rho 1998). In this paper, we present follow-up spectral 
observations to see which additional infrared lines are bright from molecular supernova shocks. These 
observations are the first detection of far-infrared H 2 O, OH, and CO lines from supernova remnants. We 
report here far-infrared spectral observations of a single position in each of 3 supernova remnants: 3C 391, 
W 44, and W 28. 


2. Observations 


All observations reported here were performed with the ISO Long-Wavelength Spectrometer (plegg et 


al. 1996| ). We used the medium-resolution grating to fully sample the wavelength range from 42.2-188.6 
fim. The coordinates are based on bright OH 1720 MHz masers ( Frail et al. 1996 ) in each remnant (for 
W 44: 18^56’"28.4'* -F01°29'59", for W 28: 18^01^52.3® -23°19'25"; B1950); for 3C 391 we shifted closer 
to the peak of shocked CS and CO broad-molecular-line emission 3C 391:BML (18^46"^47.1'* -01°00'51") 
( Reach and Rho 1998 ). The LWS beamsize is 80", and the spectra are severely affected by fringes, due 
to constructive and destructive interference (inside the spectrometer) of wavefronts from our structured, 
extended sources. We removed the fringes using the ISO Spectral Analysis Package ^ assuming that the 
emitting region is extended at all wavelengths. 


An important portion of the spectrum of 3C 391:BML, including the low-lying transition of ortho-H 20 
at 179.5 ^m, is shown in Fig. The continuum is due to dust from the supernova remnant (some 30% of 
the total) and unrelated interstellar material along the line of sight (see Paper I). In addition to the H 2 O 
and CO lines indicated, there are some remaining ripples in the spectrum that do not fall at the wavelength 
of predicted lines nor do they have the shape expected for an unresolved spectral line. A second spectral 
observation toward 3C 391:BML was performed in order to confirm some of the spectral lines, and to search 
for lines of H 2 O and OH from energy levels higher than we had already detected. The result for the OH 
(^ni_|, J = 5 “ |) line is shown in Figure]^ 


^The ISO Spectral Analysis Package (ISAP) is a joint development by the LWS and SWS Instrument Teams and Data 
Centres. Contributing institutes are CESR, IAS, IPAC, MPE, RAL and SRON. 
















































Fig. 1.— A portion of the far-infrared spectrum of 3C 391:BML obtained the ISO LWS. The combination 
of a linear baseline and 4 lines (with the instrumental line profile) is shown as a solid curve. Two H 2 O lines 
and two CO lines are labeled by species and transition. 
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3. Results 

Very bright lines from atomic and ionic C, N, and O, and weaker molecular lines were detected from 
all three supernova remnants. We will compile and present the complete line lists in a future paper. Of 
the three sources that we observed, only 3C 391:BML had a rich spectrum of molecular emission. A list 
of H 2 O, OH, and CO lines detected or limited toward 3C 391:BML is shown in Table |^. For W 28, the 
CO(16-15) and CO(15-14) emission lines were detected, and for W 44 only the CO(16-15) line was detected; 
the CO lines are a factor of 2 and 4 fainter for W 28 and W 44 compared to 3C 391:BML. The 95% 
confidence upper limit on the H20(2i2-loi) 179.5 /rm line for W 44 and W 28 is half the brightness of the 
line detected toward 3C 391:BML. Considering that 3C 391 also had the brightest [O I] 63 /rm line and 
dust continuum, it is possible that both W 28 and W 44 could have similar shock-excited spectra, with 
the H 2 O lines unfortunately just below our detection limit. While no OH emission was detected for W 28 
and W 44, the 119 /rm transition from the ground-state was detected in absorption, presumably due to 
foreground gas. One source of uncertainty is the identification of the lines. Although the resolving power 
of these observations was less than 300, such that wavelengths cannot be measured sufficiently accurately 
for conclusive identifications, our main conclusion—namely, that H 2 O, OH, and CO were detected from 
shock-excited gas in supernova remnants—is strong. In the time-honored tradition of spectroscopy, we have 
detected multiple lines from each species. 

All of the lines we believe we have detected are from the lowest accessible energy levels of each species. 
For CO, the highest level we detected was J = 16, which lies 750 K above the ground state. The energy 
level diagrams of H 2 O and OH are shown in Figures || and In both cases, we indicate the fast transitions 
whose wavelengths that fall within the range of our ISO spectrum. Including marginal detections, we see 
evidence for almost all of the fast transitions of both ortho and para H 2 O among energy levels within 320 K 
of the ground state. The H20(4i3 — 322 ) transition is masked by a bright [O I] line nearby, but there is a 
hint of the principal line from the next-highest energy level (330 — 221 ). No transitions involving levels more 
than 400 K above ground are seen. The most important emission line that is expected to be bright but is 
not detected is the 119 /rm line of OH. But this line connects directly to the ground state, and our models 
(below) predict significant foreground absorption; therefore, we suspect that the 119 ^m line is extinguished 
by foreground gas. Indeed, our spectra of W 44 and W 28, which have comparable sightlines through the 
galactic disk, show absorption features at 119 ^m (with an optical depth of order 0.1) due to foreground 
OH absorbing the dust continuum from the supernova remnants. 


4. Discussion 


4.1. Excitation and abundance of H 2 O, OH, and CO 


In order to determine the physical conditions that can produce the observed suite of spectral lines, 
and to measure the abundances of the observed species, we compared the line brightnesses to a simple 
model that balances collisional and radiative transitions within a uniform emitting region. We modeled the 
emission spectra of regions with a range of H 2 volume density and kinetic temperature, and the absorption 


spectrum of a cold slab of foreground gas with nominal molecular abundances (from Irvine et al. 1987) and 
an H 2 column density of 10^^ cm“^. We solved iteratively for the excitation, modifying radiative rates by 
the escape probability for a line profile with a width of 30 km s“^, as was found from the millimeter-wave 


CS and CO observations ( Reach and Rho 1998 ). Collision rates were taken from Offer, van Hemert, & van 
Dishoeck (1994) for OH and Green, Maluendes, & McLean (1993) for H 2 O, and the radiative transition 




Table 1. Molecular Lines from 3C 391:BML 


transition^ 

flux 

(nW m~^ sr~^) 

Ab 

(pm) 

SNR 

CO (14-13) 

9.0 ± 2.2 

186.06 

2.9 

CO (15-14)'= 

28.8 ± 2.3 

173.62 

6.4 

CO (16-15)'= 

10.0 ± 1.7 

162.86 

5.0 

CO (17-16) 

10.3 ±2.8 

153.51 

3.0 

CO (18-17)'= 

< 17 

144.78 


CO (19-18) 

< 17 

137.20 


CO (20-19) 

< 7 

130.37 


H 2 O (4i3-4o4) 

< 8 

187.11 


HfO ( 2 i 2 -loi)'=’'i 

6.6 ± 1.6 

180.69 

3.1 

H 2 O (2l2-loi) 

10.0 ±0.9 

179.59 

8.2 

H 2 O (3o3-2l2)'^ 

6.3 ±3.0 

174.54 

2.9 

H 2 O (3l3-2o2) 

< 13 

138.53 


H 2 O (423-4i4) 

< 7 

132.41 


H 2 O (4o4-3i3) 

< 13 

125.35 


H 2 O (4l4-3o3) 

< 5 

113.54 


H 20 ( 221 - 110 ) 

< 7 

108.07 


H 2 O (220-1ii)‘= 

7.1 ± 2.0 

100.47 

2.9 

H 2 O (5l5-4o4) 

< 6 

95.63 


H 2 O (322-2ii) 

< 7 

89.99 


H 2 O (321-2i2) 

< 8 

75.38 


OH (1-1;3-1)'= 

4.2 ± 1.2 

163.65 

2.6 

OH (3-3;5-3) 

< 7 

119.33 


OH (l-3;5-7) 

< 13 

115.22 


OH (l-l;5-3) 

3.6 ± 1.4 

98.69 

2.6 

OH (l-3;3-5) 

< 9 

96.34 


OH p-3;7-5) 

12.7 ±3.3 

84.58 

3.1 

OH (l-3;l-3) 

17.4 ± 1.9 

79.11 

11.9 

OH (3-3;9-7) 

< 13 

65.20 


OH (l-3;3-3) 

< 13 

53.30 



^CO transitions are labeled in the format (a-6) where 
a and b are the upper and lower rotational quantum 
number; OH transitions are labeled in the format {a-b;c- 
d) indicating the transition from level ^ to 

level J — d/2. 

^observed line center for detected lines, accurate to 0.2 
(0.1) ^m for A greater (less) than 88 pm; laboratory rest 
wavelength for undetected lines 

^line potentially blended with another 

"^This line could be the H 2 O ( 22 i- 2 i 2 ) transition 
instead, but the excitation model suggests this transition 
contributes less than 10% of the observed brightness. 
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rates were taken from Pickett et al. (1996). The observed brightness ratio of 79 to 84 /rm lines of OH is 
sensitive to the gas density, suggesting n(H 2 ) ~ (1 — 6) x 10® cm“®. The ratios of the 5 ^ 4, 3 —> 2, and 
2 —> 1 millimeter lines of CS are also sensitive to the gas density, suggesting n(H 2 ) ~ (3 — 4) x 10® cm“®, 
and the millimeter CS and CO lines somewhat constrain the temperature, T > 50 K ( Reach and Rho 1998 ). 
For the abundance calculations, we will assume n(H 2 ) = 2 x 10® cm ® and 100 < T < 1000 K, which is 
consistent with the presence and lack of other OH and H 2 O lines in the observed wavelength range. 


(Reach and Rho 1998 


We determine the abundances of the various molecules assuming all lines arise from the same physical 
region, with constant temperature and density. The column density of OH is ~ 2 x 10^® cm“^, and the 
optical depth of the 79 /rm line is of order unity. The CO excitation is much more sensitive to temperature; 
to match the brightness of the well-detected 15 ^ 14 line, the CO column density ~ 2 x 10^®(r/200)“® 
cm“^. The H 2 O lines are estimated to be optically thick: if we associate the line observed at 180.69 p,m with 
H2®0 (2 i2 -1 io), then the isotope ratio implies an optical depth ^ 400 for the brightest H 2 O line. However, 
in the low-density limit, n(H 2 ) < 10® cm“®, spontaneous decay is still faster than collisional de-excitation, 
and the line brightness still measures the H 2 O abundance ( Irvine et al. 1987 ). The column density of H 2 O 
from the excitation model is ~ 3 x 10^^ cm“^. To determine the absolute abundance of each species, we 
require a measure of the H 2 column density. If we assume the emitting region is a uniform sphere with 
diameter equal to the observed angular size of shock-excited molecular gas (^ 

), and we take from the excitation model n(H 2 ) = 2 x 10® cm“® 


30" from [Reach and Rho 
we find 7V(H2) 


8 X 10^3 cm-3. This 

value of 1V(H2) agrees with the observed brightness (from our ISO SWS observations, in preparation) of the 
S(3) line of H 2 , if T ~ 200. In summary, we estimate the relative abundances of C0:0H:H20 to be 100:1:7, 
and the abundance of water in the shocked cloud is [H20/H2]~ 4 x 10“^. 


4.2. Chemistry of H 2 O and OH 


The warm H 2 O in 3C 391:BML is likely due to shock-enhanced chemistry. In relatively hot gas, OH 
is rapidly converted into H 2 O by the reaction OH -|- H 2 ^ H 2 O -I- H, so that all the available O that is 
not already locked in CO would be converted into H 2 O (Draine et al. 1983, McKee & Hollenbach 1979). 
Models for the oxygen chemistry predict very efficient conversion of OH into H 2 O for non-dissociative 
shocks into high-density (ng = 10® cm”®) clouds but comparable OH and H 2 O column densities for shocks 
into intermediate-density {no = 10® cm”®) clouds, while CO and H 2 O are comparably abundant (Graff & 
Dalgarno 1987 ). Behind a fully dissociative shock, the molecules re-form in a much cooler (~ 100-500 K) 
region ( Hollenbach fc McKee lO"^ ). In cooler gas, the 1420 K barrier ( Wagner fc Graff 1987 ) for the 
OH -I- H 2 ^ H 2 O -I- H reaction cannot be overcome, and OH is more abundant than H 2 O (McKee & 


Hollenbach 1979, van Dishoeck and Black 1986) except deep inside dense molecular cores ( ^ternberg &j 
Dalgarno 1995). The chemistry in cooler gas depends on photo-dissociation; OH is much more abundant 
than H 2 O ( van Dishoeck and Black 1986| ), except deep inside very dark cores, where the photodissociation 
rate is very low (Bergin et al. 1995). 


The relative strengths of the H 2 O lines that we observe were calculated by Kaufman & Neufeld (1996), 
and our observations are generally consistent with the models for pre-shock density ng = 10^”® cm”®. 
However, these models predict a very low OH abundance: the H 2 O 179.5 /rm line is predicted to be two 
orders of magnitude brighter than the OH 84.5 /rm doublet, while we observe comparable brightnesses. The 
relatively low excitation we observe for the H 2 O and OH molecules {Tupper < 400 K) suggests that the gas 
we are observing is not presently hot enough for the rapid H 2 O production. But this does not preclude 
the H 2 O having formed in a short-lived, high-temperature region, and we observe only the cooling region. 






































Fig. 3.— Energy-level diagrams for ortho (left) and para (right) H 2 O. All levels within 500 K of the 
ground state are shown. The fastest radiative de-excitations of each level are indicated with arrows, if 
the wavelength of the transition is within the range of our ISO spectrum. The arrows were drawn with 
the following characteristics: detected lines have a solid-filled arrowhead, possible detections have an open 
arrowhead, and non-detections have dashed arrow. 
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Fig. 4.— Energy-level diagram of OH. Each level is labeled by its rotational quantum number J. For the 
transitions within the range of our ISO spectrum, arrow symbols are the same as in Fig. Clearly detected 
lines have a solid-filled arrowhead, possible detections have an open arrowhead, and non-detections have 
dashed arrow. The transitions are labeled with the wavelength in /im. 
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We observe more H 2 O than OH, suggesting that high-temperature chemistry was operative long enough to 
leave a lasting effect on the chemistry of this gas. 


The detection of OH, with an abundance only 7 times less than that of H 2 O, disagrees with theoretical 
models for C shocks, which predict nearly complete conversion of OH into H 2 O. This suggests either (1) the 
H 2 O abundance is underestimated due to beam-dilution, (2) we are observing OH and H 2 O from different 
types of shocks in regions with a range of pre-shock densities, or (3) the models are not appropriate for 
the molecular shocks in 3C 391. We detected bright ionic lines (including [O HI] and [N HI]) that are 
only expected from dissociative shocks, so we know that a range of shocks is present within our beam. 
However, the infrared OH line ratios are consistent with coming from a region with the same physical 
conditions as inferred from the infrared H 2 O, millimeter-wave CS, and radio-wave OH maser lines. The 
disagreement between observations and predictions is at least partially due to the fact that the models are 


for non-dissociative shocks (Kaufman and Neufeld 1996). Including the effects of photodissociation of H 2 O 
by the interstellar radiation field—as well as radiation local to the remnant— could produce OH in the 
abundance we observe (Lockett, Gauthier, & Elitzur 1998). 

The water abundance we derive for 3C 391:BML is significantly lower than that observed from 
molecular shocks in HH 54 ( Liseau et al. 1996 ) or Orion ( |Harwit et al. 19^ ). In particular, in the Orion 
BN-KL region, the H 2 O abundance is [H 2 0/H2]ci 5 x 10“^ ( Harwit et al. 1998 ), incorporating nearly all of 
the oxygen in the gas. The difference between the Orion BN-KL shock and the 30 391:BML shock could 
be due to limited angular resolution: 30 391 is about 20 times further away than Orion, while the angular 
resolution of our observations is the same as that of Harwit et al. (1998). The difference could also be 
due to details of the interaction between the energetic events (a steady stellar wind in Orion vs. impulsive 
supernova shock in 3G 391) and the interstellar medium (a molecular cloud and H H region in Orion vs. a 
giant molecular cloud with no H H region in 3G 391). 


5. Conclusions 


For the first time, thermal emission from the lower energy-levels of H 2 O and OH were detected from 
a supernova remnant. The emission arises from gas cooling behind a shock front that is impinging a 
particularly dense clump in the parent molecular cloud; this site is called 30 391:BML (broad molecular 
line). Our spectra of interaction sites in W 44 and W 28 did not reveal OH or H 2 O emission, although 
CO emission and OH (foreground) absorption were detected. 30 391 is not a unique case of a supernova 
remnant-molecular cloud interaction, but the infrared molecular lines imply the presence of higher-density 
clumps than near W 44 and W 28, similar to 10 443. Together with outflows from young stars, supernova 
remnants are energetic events capable of altering cloud chemistry and producing substantial columns of H 2 O 
and OH. Supernova-molecular cloud interactions are probably common, and they have been suggested as an 


explanation for a new class of mixed-morphology X-ray and radio supernova remnants ( Rho fc Petre 1998 ). 
We can expect substantial advances in understanding shock-induced chemistry in supernova-molecular 
cloud interactions by using future observatories capable of higher angular and spectral resolution in the 
far-infrared, such as planned for SOFIA ( Erickson 1995| ) and FIRST (Poglitsch 1998), where we will better 
resolve the fast and slow shocks into gas of varying pre-shock density, and we will better resolve a larger 
number of spectral lines from each molecule. 


The observations that we describe in this paper are part of the ISO open time granted to U.S. 
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